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In our recent Letter [l[ , we reported a study of how the 
thermal conductivity of YBa2Cu30j,(YBCO) changes as 
its doping is made to cross the superconducting critical 
point at p c = 0.05, i.e., from pi < p c to pf > p c . For the 
same crystal with the same contacts, we showed that the 
difference Ak = n(pi) — n(pf) = (3T 3 up to 500 mK. This 
is an experimental fact, free from any analysis and con- 
firmed on several specimens. We attribute this difference 
to heat-carrying bosons with a T 3 conductivity analogous 
to the magnons in the undoped cupratc antiferromagnet 
Nd 2 Cu0 4 0). 

Twin-boundary scattering. - Ando [3J points out cor- 
rectly that YBCO undergoes a structural transition from 
a tetragonal (tetra) phase for p < p c to an orthorhom- 
bic (ortho) phase for p > p c . He raises the possibil- 
ity that in twinned samples like ours, scattering from 
twin boundaries will cause a drop in phonon conductiv- 
ity when crossing into the ortho phase. While this may 
be true, it does not explain the [3T 3 difference we observe 
because phonon conductivity in (twinned or untwinncd) 
cuprate crystals is never T 3 , at any doping. This was 
shown very clearly in our study of Nd2CuC>4 Q, where 
the magnon conductivity goes precisely as T 3 , while that 
of phonons goes as T 26 . Below 0.5 K or so, the mean 
free path of magnons is limited by the rough sample edges 
(quasi-2D) while that of phonons is limited by the sam- 
ple faces (3D) which, in as-grown crystals, are mirror-like 
and cause specular reflection. This introduces a strong In- 
dependence to the phonon mean free path, which means 
that K, p is no longer cubic in T, as shown in our recent 
study of as- grown and roughened samples of Nd2CuC>4 
Q. Therefore, the T 3 term we observe in YBCO cannot 
be due to a change in n p caused by additional phonon 
scattering for p > p c . If significant at all, this additional 
scattering will be seen at higher temperature, above the 
boundary scattering regime (T > 0.5 K or so). (Note 
that in their own measurements on YBCO, Ando and 
co-workers claim to see a phonon T 3 term at low T 
However, this is always limited to temperatures below 
150 mK or so.) 

Although their data is not shown in our Letter 
we have measured samples whose entire doping evolution 
took place in the tetra phase, i.e., below 5%. We refer to 
these samples (labeled M, N, O) in our Letter and in FigQ] 
we show data for the most underdoped one (sample O, 
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FIG. 1: Top: Boson - Change in the thermal conductivity 
k of YBCO sample O (y = 6.31) between po = 4.3% and 
pi = 4.6%, well below the I-S transition at p — 5%. The 
reduction in k, giving Ak oc T 3 , occurs entirely within the 
tetragonal phase and cannot be due to scattering off twin 
boundaries. Inset: sketch of how k would evolve if twin 
boundary scattering of phonons were causing a reduction 
(dashed line) and how k is observed to change in our mea- 
surements (solid line) (see Fig. 2 in [l[). Bottom: - Change 
in k of YBCO sample L (y = 6.33) between po — 4.7% and 
pi — 5.0%, with the (3T 3 contribution subtracted, plotted as 
Ak/T — f3T 2 vs T. The corresponding change in electrical 
conductivity a (see Fig. 3 of [J|) is plotted as LoAa vs T. 
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with y = 6.31). The change in conductivity between po = 
4.3% and pi = 4.6% exhibits the same T 3 dependence 
shown in our Letter [l[ . Because sample O never entered 
the ortho phase, twin boundaries cannot be invoked as 
the cause of the reduction in k. Invoking strains that 
build up as the tetra-ortho transition is approached from 
below also does not work, as scattering from strain fields 
is known to vary linearly with phonon frequency and so 
causes k p oc T 2 , not T 3 . 

Wiedemann- Franz law. - In our Letter, we fit the mea- 
sured k(T) to a sum of three terms: k(T) = aT + /3T 3 + 
r(T) coming from fermions, bosons, and phonons, re- 
spectively. Ando argues that we should have explicitly 
included an additional term related to charge conduction. 
However, the only way one can do this a priori is to use 
the Wiedemann-Franz (WF) law: n e — L Ta = L T/p, 
where Lo is the Sommcrfcld constant and a the electrical 
conductivity. In other words, Ando wants us to replace 
Kf = aT by Kf — a'T + LqT/ p. 

Moving away from arguments over fitting procedure, 
the best way to compare charge and heat conductivities 
is to directly examine the respective changes induced by 



doping. In the bottom panel of Figfl] we plot the differ- 
ence in k for sample L between po = 4.7% and p\ = 5.0%, 
without the T 3 term, i.e., Ak/T - (3T 2 vs T. We com- 
pare this with the corresponding change (on the same 
sample at the very same dopings) in electrical conductiv- 
ity, L Q Aa [l[. Because k/T does not change with doping 
while a does, the WF law is unambiguously violated, in- 
dependently of the fitting precedure. This shows that 
Ando's suggestion of including an explicit n e term leads 
to an inconsistency, in the sense that it relies on the WF 
law which does not hold here. 



* Electronic address: Louis. Taillefer@UShcrbrooke.ca 
[1] N. Doiron-Leyraud et al, Phys. Rev. Lett. 97, 207001 
(2006). 

[2] S.Y. Li et al, Phy s. Rev. Lett. 95, 156603 (2005). 
[3] Y. Ando, preprint arXiv:0711.4654vl (2007). 
[4] S.Y. Li et al., preprint arXiv:0709.3075 (2007). 
[5] X.F. Sun et al, Phys. Rev. Lett. 93, 107001 (2004). 



